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FOURTH QUARTER SECOND YEAR PROGRESS REPORT 

This report covers technical progress during the fourth quarter of the second year of 
NASA Sun-Earth Connections Theory Program (SECTP) contract “The Structure and 
Dynamics of the Solar Corona and Inner Heliosphere,” NAS5-99188, between NASA and 
Science Applications International Corporation, and covers the period May 16, 2001 to August 
15, 2001. Under this contract SAIC and the University of California, Irvine (UCI) have 
conducted research into theoretical modeling of active regions, the solar corona, and the inner 
heliosphere, using the MHD model. 

Press Release at the Spring AGU Meeting, Boston, MA, May 30, 2001 

Zoran Mikic was invited by the AAS/SPD press officer. Dr. Criag DeForest, to give a 
presentation at a press conference about the total solar eclipse that was expected in Africa on 
June 21, 2001. An accompanying press release is contained in Appendix A. 

Presentations at the Spring AGU Meeting, Boston, MA, May 29— June 2, 2001 

Our progress during this reporting period is summarized by the following papers that 
were presented by SAIC staff members at the Spring AGU Meeting, Boston, MA, May 29- 
June 2, 2001. 

Predicting the Structure of the Solar Corona During the 21 June 2001 Total Solar Eclipse 
Z. Mikic, J. A. Linker, R. Lionello, and P. Riley (SAIC, San Diego) 

We describe the application of a three-dimensional magnetohydrodynamic (MHD) model 
to the prediction of the structure of the corona during the total solar eclipse that is expected to 
occur on 21 June 2001. The calculation uses the observed photospheric radial magnetic field as 
a boundary condition. This model makes it possible to determine the large-scale structure of 
the magnetic field in the corona, as well as the distribution of the solar wind velocity , plasma 
density, and temperature. We will use magnetic fields observed on the solar disk prior to 
eclipse day to predict what the corona will look like during the eclipse. The estimated coronal 
density will be used to predict the plane-of-sky polarization brightness prior to the eclipse. 

A copy of this presentation appears in Appendix B. 


Are There Two Classes of Coronal Mass Ejections? A Theoretical Perspective 
J. A Linker, Z. Mikic, R. Lionello, and P. Riley (SAIC, San Diego) 

Coronal mass ejections (CMEs) are generally accepted as the cause of nonrecurrent 
gnetic storms at Earth. Statistical compilations of CME events have shown that CMEs 
launched in the corona can have a wide variation in speeds (Hundhausen et ai, JGR 99, 6543, 
1994). The speed of the CME at Earth and the presence or absence of an interplanetary 
shockwave is an important component of the geoeffectiveness; therefore, the mechanism(s) by 
which fast CMEs might be produced are considered to be of particular importance. Recently, 
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the examination of the acceleration profiles of CMEs has led to the possible classification of 
CMEs as either constant speed CMEs or constant acceleration (St. Cyr et ai, JGR 104, 12493, 

1999; Sheeley et al., JGR 104, 24739, 1999). In this talk, we will examine reasons why one 
might expect CMEs to show two classes of acceleration profiles, and if they do, what the 
implications are for models of CME initiation. Examples from MHD simulations of CMEs for 
different initiation mechanisms will be used to illustrate essential points. 

^Research supported by NASA and Boston University’s Integrated Space Weather Modeling 
project (funded by NSF). 

A copy of this presentation appears in Appendix C. 


Using Global MHD Simulations to Interpret In Situ Observations of CMEs 

P. Riley, J. A. Linker, R. Lionello, Z. Mikic (SAIC, San Diego) 

D. Odstrcil, V. J. Pizzo (SEC, Boulder) 

T. H. Zurbuchen (U. Michigan) 

D. Lario (JHU/APL) 

In this study, we combine two MHD models to simulate the initiation, propagation, and 
dynamic evolution of flux-rope-like CMEs through the corona and out to 1 AU. The coronal 
model encompasses the region of the solar corona from 1 R s to 20 R s , while the heliospheric 
model encompasses 20 R s to 1 AU. The CME initiated in the corona propagates smoothly 
across the outer boundary of the coronal solution and through the inner boundary of the 
heliospheric solution. The model solutions show a rich complexity, which, given the relative 
simplicity and idealization of the input conditions, bear a strong resemblance to many observed 
events, and we use the simulation results to infer the global structure of some of these 
observations. In particular, we highlight an event that was observed by both Ulysses and ACE 
in February/March, 1999. At this time, Ulysses was located at - 5 AU and 22°S heliographic 
latitude; thus the two spacecraft were separated significantly both in heliocentric distance and 
latitude. We also use these simulations to separate dynamical effects from force-free models of 
flux ropes in the solar wind. 

A copy of this presentation appears in Appendix D. 


Modeling ofTransequatorial Loops with MH4D 
R. Lionello and D. Schnack (SAIC, San Diego) 

MH4D (Magnetohydrodynamics on a TETRAhedral Domain) is a new algorithm to 
perform magnetohydrodynamic (MHD) simulations of the active regions of the Sun, including 
the large scale coronal structure that surround them. MH4D is a massively-parallel, device- 
independent numerical code for the advancement of the resistive and viscous MHD equations 
on an unstructured grid of tetrahedra. The use of an unstructured grid allows us to increase the 
resolution in the regions of physical interest. A variational formulation of the differential 
operators ensures accuracy and the preservation of the analytical properties of the operators 
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(V-B = 0, self-adjointness of the resistive and viscous operators). The combined semi-implicit 
treatment of the waves and implicit formulation of the diffusive operators can accommodate the 
wide range of time scales present in the solar corona. The capability of mesh refinement and 
coarsening is also included. A preliminary result is presented: a simulation of transequatorial 
loops that include fine details of two interconnected active regions. 

A copy of this presentation appears in Appendix E. 
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APPENDIX A 


Press Release at the Spring AGU Meeting 
Boston, MA, May 30, 2001 
Presented by Zoran Mikic 
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PREDICTING THE STRUCTURE OF THE SOLAR CORONA 
DURING THE 21 JUNE 2001 TOTAL SOLAR ECLIPSE 

Total solar eclipses offer an excellent opportunity to observe the solar corona. During a total solar 
eclipse the moon blocks the bright light from the solar disk, allowing the faint light scattered by the solar 
corona, which is more than a million times fainter than the photosphere, to become visible. During totality 
the structures that characterize the white-light corona become apparent, including prominences, helmet 
streamers, polar plumes, and coronal holes. Observers who witness a total solar eclipse invariably report 
that it is a beautiful sight to behold. 

On 21 June, 2001 a total eclipse of the Sun will be visible in the southern hemisphere, beginning in 
the South Atlantic, crossing southern Africa and Madagascar, and terminating in the Indian Ocean. Drs. 
Zoran Mikic, Jon Linker, Pete Riley, and Roberto Lionello, of Science Applications International 
Corporation (SAIC) in San Diego, California, have developed a theoretical model to predict what the solar 
corona will look like during forthcoming total solar eclipses. Their model has been used to predict the 
shape of the corona during the eclipse that is expected to occur on 21 June. 

The model is based on the three-dimensional magnetohydrodynamic (MHD) equations that describe 
the interaction of the solar wind with coronal magnetic fields. The group’s results, which have been 
financially supported by the National Aeronautics and Space Administration (NASA) and the National 
Science Foundation (NSF), are being presented by Dr. Mikic at the American Geophysical Union meeting 
in Boston. 

The calculation relies on Earth-based measurements of the magnetic field in the solar photosphere to 
infer the structure of the solar corona. The measurements are taken at the National Solar Observatory at 
Kitt Peak . “It is remarkable that measurements of the magnetic field in the photosphere can tell us so 
much about the corona,” says Dr. Mikic. “The simulations already have a strong resemblance to coronal 
images. In the future, the agreement will only improve as we refine the physics in our model and as we 
utilize ever-faster computers. We will be able to study even finer details in the corona.” The researchers 
use supercomputers at the San Diego Supercomputer Center (part of NSF’s National Partnership for 
Advanced Computing Infrastructure) and the Department of Energy’s National Energy Research 
Supercomputer Center to solve the equations. The output from the model is used to predict the brightness 
of the corona. 

So far, this model has been applied to five eclipses (see http://haven.saic.com) . This kind of 
modeling helps us to understand the structure of the solar corona, especially the location of helmet 
streamers and coronal holes, and the nature of the fast and slow solar wind, and provides a rudimentary 
test of predictive capability. Eclipses offer an opportunity to test such models, and to understand the 
influence of the Sun on the Earth. 



Predicted polarization brightness (pB) for the 21 June, 2001 eclipse, together with traces of the magnetic field lines, at 13: 10 
UT (corresponding to totality in Lusaka, Zambia). The pB signal is produced by white light scattered off electrons in the 
coronal plasma. This is the view of the Sun that would be seen by an observer on Earth with a camera aligned so that 
vertical is toward the Earth’s north pole. 

Zoran Mikic, SAIC, 10260 Campus Point Drive, San Diego, CA 92121; Tel: 858-826-6934; Email: mikicz@saic.com 
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APPENDIX B 


“Predicting the Structure of the Solar Corona During the 21 June 2001 Total Solar Eclipse” 
Z. Mikic, J. A. Linker, R. Lionello, and P. Riley 
Presented at the Spring AGU Meeting 
Boston, MA, May 29-June 2, 2001 
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Equations The Polytropic Model 

| • Neglect thermal conduction, coronal 
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MHD Equations Finding MHD Solutions 

(Polytropic Model) TI c . , 

• Use line-of-sight magnetograms to 
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Polarization Brightness ! "Canonical" Helmet Streamer 

Light scattered off the coronal electrons • Idealized helmet streamer configuration 
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August 11, 1999 eclipse (Central 
Europe, Turkey, Iran): 

We predicted the coronal structu 
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Predicted on May 19, 2001 using magnetic field data from Kilt Peak Solar 
Observatory for Carrington rotation 1975 (April 9 - May 6, 2001). 
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APPENDIX C 


“Are There Two Classes of Coronal Mass Ejections? A Theoretical Perspective” 
J. A Linker, Z. Mikic, R. Lionello, and P. Riley 
Presented at the Spring AGU Meeting 
Boston, MA, May 29-June 2, 2001 
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My Assignment:" Questions 

The modelers should address the question • Does the division of CMEs into 
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Eruption of a 3D Helmet Streamer triggered Shearing Flows 



Consider 5 Model Problems: 
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Relaxed Helmet Streamer 





Eruption of a Helmet Streamer Eruption of a Helmet Streamer 

By Emerging Flux By Emerging Flux 

Flux Y(r,z) Polarization Brightness 
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TimeHeight Profiles: Flux Cancellation and Shearing Simulations 


Flux Cancellation 1 


Shearing Motions 1 



Time (Hours) 


Time (Hours) 




Flux Cancellation, Medium Wind Background 
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What Have We Learned? 


q j 

CD 

03 

u 


c o 

d 

o 


00 "73 

4 -> 

3 

a 


.3^ O 

S Cl, £ 

C/3 CO 


u 

M-| 

o 

d 

o 


d g o 

£ 
g 


CD . 

■i-H "» 


.£'2 
'V g 
03 

° 

O.S 

5-1 4-> 


a; 

d 

o 

d f d 

03 d 


V 2 

r-H QJ 


‘43 

CD <3 
M *j3 

d'2 


to 'pi^TtS 

s y 

q^Io 


CD 

n3 

CD 

w 


u 


c^* 

S 

CD 

jS ^ 

<D 

*0,5 
a a 

d 

d o 
.2^ 
.2^ 

> o 

CO^ 

;d73 

s & 

CO 03 

gS 6 

Q.SZ 



8 w-6 

gs g 

co fj £ 
co ^ ~ 
03 73 ^ 
u d xi 
co - 2 ’co 

•d CO CO 

rS ^ O 
ti d< 

* d_2 

U 3-5 



r 2 So w 

ss 

8 

a.ft’d -a 
“ CT3 g 
■g g d >3 

4 -> CT* (/J 

2 ^ ±3 d3 
d co 33 73 

O d Cl, j_, 

u o d o 

5 U -d -4H 


P4 

< 

2 

2 

p 

CD 


CO 


CO 


5 d o 
’•3 Ocd 


73 O 
”d QJ co 
0) QJ Qj 

73 0 -,x: 
-d co tJ 
.» oi 
Xd dn 

CO JS 

W 2 

S§ 


o 

5-1 

d, 

0> 


U cj^d-d 
^d 
o d - d 
o £ 


d 

o 


03 


co 


W3 

.S 

IS 5 

u sn'd 

d aj d . 

U gj dS 

x u Oi *d 

3 *73-g 

* £2 
^pd U 

gx.vj* 

H cu 57*43 


os 


u -d 
•£ 

QJ d cO 
-dd 52 
Cl, g 

Jjj u 3 


o 

CD 

i—H 


.23> 

d 2 


« -2 
d u 
o (U 

X73 co *J3 ‘5T 

77 ^ 73 

73 S o 3 o> 

0) a 33 00 2^ 

4 -< M_| . 7 * rv 

.2 >^bC?B co 
43 co a S 4 -. 

§ >•907 
i> M U d 
O 03 _ < « 

co 33 a 2 "co 

5 ffw-g § 

f\ .5 . -d u 
CJ J2 M 73 
a> g u 
rdpS-Si^-Si 
H a w o >> 



SAIC-01/801 3:APPAT-280 
August 15, 2001 


APPENDIX D 

“Using Global MHD Simulations to Interpret In Situ Observations of CMEs” 
P. Riley, J. A. Linker, R. Lionello, Z. Mikic, D. Odstrcil, V. J. Pizzo, 

T. H. Zurbuchen, and D. Lario 
Presented at the Spring AGU Meeting 
Boston, MA, May 29-June 2, 2001 



Using Global MHD Simulations to 
Interpret in situ Observations of CMEs 
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Overview 


■ Introduction 

• Coronal observations of CMEs 

• In situ observations of magnetic clouds 

• Modeling Flux ropes 

■ MHD modelling of CMEs 

• Coronal model (SAIC) 

• Heliospheric model (NOAA/SEC) 

■ Comparison with observations 

• ACE/Ulysses spring 1999 event 

■ Summary 

■ Future work 




Introduction: White Light Observations 

of CMEs 



Introduction: In situ observations of 

magnetic clouds 



Introduction: Models of Magnetic 
clouds/flux ropes 


• Kinematic Models: 

■ Force-free 

* Force-free with expansion 

* None force-free with expansion 

* Dynamic Models: 

* Fluid 
- MHD 


MHD Simulation: Coronal solution 



*' <R S ) 


MHD Simulation: Coronal solution 



MHD Simulation: Coronal solution 



MHD Simulation: Heliospheric solution 



MHD Simulation: Heliospheric solution 






MHD Simulation: Heliospheric solution 



Comparison of an observed magnetic cloud 
with the generic simulation 




( Adapted from Marubashi. 1999) 








Ulysses/ACE May 1999 CME 
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Force-Free Fitting at ACE and Ulysses 


Parameter 

ACE 

Ulysses 

Start (day hr:min) 

49 14:00 

62 20:00 

Stop (day hr:min) 

50 11:00 

64 22:00 

<l> (deg) 

282.1 

271.1 

9 (deg) 

-1.3 

53 

v/R 

0.738 

0.0064 

Helicity 

-1 

-1 

R{ AU) 

0.28 

0.34 


ACE/Ulysses Feb 1999 CME Event: 
Global structure from force-free fitting 



Ulysses 









Comparison of in situ observations 
with MHD simulation: Ulysses 



Comparison of in situ observations 
with MHD simulation: ACE 





Deceleration of CME 


■ Estimate time for CME to travel from ACE to Ulysses 
Ecme(ACE) = 590 km/s 

E CM E(Ulysses) = 460 km/s 

Assume constant deceleration between 1 and 5 AU 
=> At = 13.62 days 

■ CME actually arrived 13.9 days later 

=> Large-scale structure not significantly distorted 


ACE/Ulysses Feb 1999 CME Event: 
Global structure 




Summary 


■ Idealized MHD simulations: 

• provide global context for interpreting in situ CME 
observations 

• Emphasize significant role of dynamics on the 
evolution of ejecta and their associated disturbances 

• Identify new effects that can be searched for in 
observations 


Future Studies 


■ Explore effects of different launch profiles 

• Do different mechanisms of eruptions lead to different 
interplanetary signatures? 

■ Launch CME in three dimensions 

• 3-D flow 1C 

• Magnetic connectivity to Sun 

■ Incorporate composition into model 
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APPENDIX E 


“Modeling of Transequatorial Loops with MH4D” 
R. Lionello and D. Schnack 
Presented at the Spring AGU Meeting 
Boston, MA, May 29-June 2, 2001 
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FINITE VOLUME METHOD THE MAGNETIC FIELD 
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DIVERGENCE OF B ALTERNATE DERIVATION OF B 
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THE CURRENT DENSITY SELF-ADJOINTNESS 
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VARIATIONAL PRINCIPLE |j DISCRETE MINIMIZATION 

• Minimize /(A) on a grid with N v vertices and N x 
Minimize functional tetrahedra 
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